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We investigate a 3+1 dimensional toy model that exhibits 
spontaneous breakdown of chiral symmetry, both in a light- 
front (LF) Hamiltonian and in a Euclidean Schwinger-Dyson 
(SD) formulation. We show that both formulations are com- 
pletely equivalent — provided the renormalization is properly 
done. For the model considered, this means that if one uses 
the same transverse momentum cutoff on the SD and LF for- 
mulations then the vertex mass in the LF calculation must 
be taken to be the same as the current quark mass in the SD 
calculation. The kinetic mass term in the LF calculation is 
renormalized non-trivially, which is eventually responsible for 
the mass generation of the physical fermion of the model. 



I. INTRODUCTION 

Light-front (LF) coordinates are natural coordinates 
for describing scattering processes that involve large mo- 
mentum transfers — particularly deep inelastic scatter- 
ing. This is because correlation functions at high ener- 
gies are often dominated by the free quark singularities 
which are along light-like direction. This is one of the 
main motivations for formulating field theories in the LF 
framework 0] . But light-front field theories have another 
peculiar feature, namely naive reasoning suggests that 
the vacuum of all LF Hamiltonians is equal to the Fock 
vacuum ^J,[l] . It thus appears as if LF Hamiltonians (at 
least those without the so-called zero-modes, i.e. modes 
with k + = 0) cannot be able to describe theories where 
the vacuum has any nontrivial features, such as QCD — 
where chiral symmetry is believed to be spontaneously 
broken. Even if one is only interested in parton distri- 
butions, one might be worried about using a framework 
where the vacuum is just empty space to describe a the- 
ory like QCD. 

However, it is not quite so easy to dismiss LF field the- 
ory as the following few examples show: One of the first 
field theories that was completely solved in the LF formu- 
lation was QCDi + i(Nc — » oo) ||. 't Hooft's solution did 
not include any zero-modes and therefore he had a triv- 
ial vacuum. Nevertheless, his spectrum agreed perfectly 
well with the numerical results from calculations based 
on equal time Hamiltonians [0. Beyond that, applica- 
tion of current algebra sum rules to spectrum and decay 
constants obtained from the LF calculation formally gave 
nonzero values for the quark condensates that also agreed 
with numerical results at equal time pj. This peculiar 



result could be understood by defining LF field theory 
through a limiting procedure, which showed that some 
observables (here: spectrum and decay constants) have a 
smooth and continuous LF limit, while others (here quark 
condensates) have a discontinuous LF limit. Other ex- 
amples have been studied, in which it was still possible 
to demonstrate equivalence between LF results and equal 
time results nonperturbatively, provided the LF Hamilto- 
nian was appropriately renormalized . Even though 
these examples are just 1+1 dimensional field theories, 
it is generally believed among the optimists in the field 
[jlT3| , P~l| Jlj that it should be possible in 3+1 dimensional 
field theories as well to achieve equivalence between spec- 
tra of LF Hamiltonian and equal time Hamiltonians by 
appropriate renormalization. However, no nontrivial ex- 
amples (examples that go beyond mean field calculations) 
to support such a belief existed so far. 

In this paper, we will give a 3+1 dimensional toy model 
that can be solved both in a conventional framework 
(here by solving the Schwinger-Dyson equations) as well 
as in the LF framework. We will unashamedly omit zero 
modes as explicit degrees of freedom throughout the cal- 
culation. Nevertheless, we are able to show that appro- 
priate counter-terms to the LF Hamiltonian are sufficient 
to demonstrate equivalence of the spectrum and other 
physical properties between the two frameworks. 



II. A SIMPLE TOY MODEL 

The model that we are going to investigate consists of 
fermions with some "color" degrees of freedom (funda- 
mental representation) coupled to the transverse compo- 
nent of a vector field, which also carries "color" (adjoint 
representation). The vector field does not self- interact .[] 
Furthermore, we will focus on the limit of an infinite 
number of colors: Nc — > oo (g fixed), which will render 
the model solvable in the Schwinger-Dyson approach 



1 Note that, for finite Nc, box diagrams with fermions would 
induce four boson counter-terms, which we will ignore here 
since we will consider the model only in the large Nc limit. 
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With "_L component" we mean here the x and y compo- 
nents. Furthermore we will impose a transverse momen- 
tum cutoff on the fields and we will consider the model 
at fixed cutoff. Also, even though we are interested in 
the chiral limit of this model, we will keep a finite quark 
mass since the LF formulation has notorious difficulties 
in the strict m = case. Those difficulties can be avoided 
if one takes m > considers m — > 0. 

Even though certain elements of the model bear some 
resemblance to terms that appear in the QCD La- 
grangian, the model seems is a rather bizarre construc- 
tion. However, there is a reason for this: What we are 
interested in is a LF-investigating of a model that ex- 
hibits spontaneous breakdown of chiral symmetry. Fur- 
thermore, we wanted to be able to perform a "refer- 
ence calculation" in a conventional (non-LF) framework. 
In the large Nc limit, the rainbow approximation for 
the fcrmion self-energy becomes exact, which allows us 
to solve the model exactly in the Schwinger-Dyson ap- 
proach. The vector coupling of the bosons to the fermions 
was chosen because it is chirally invariant and because a 
similar coupling occurs in regular QCD. The restriction 
to the _L component of the fields avoids interactions in- 
volving "bad currents" . Finally, using a transverse mo- 
mentum cutoff both in the Schwinger-Dyson approach 
and in the LF calculation should allow us to directly com- 
pare the two formulations. 



A. Schwinger-Dyson Solution 

Because the above toy model lacks full covariance 
(there is no symmetry relating longitudinal and trans- 
verse coordinates) the full fermion propagator is of the 
form 

S F ( P n =l/ L S L (f L ,f ± )+ fi±S ± (? L ,? ± ) + S (f L ,f ± ), 

(2.2) 

where }/l = fco7° + ^37 3 and f/± = fci7 1 + &27 2 - On very 
general grounds, it should always be possible to write 
down a spectral representation for Spf\ 



2 What we need is that the Green's functions are analytic ex- 
cept for poles and that the location of the poles are consistent 
with longitudinal boost invariance (which is manifest in our 
model) . The fact that the model is not invariant under trans- 
formations which mix pl and p± does not prevent us from 
writing down a spectral representation for the dependence on 
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where i — L, _L, 0. Note that this spectral representation 
differs from what one usually writes down as a spectral 
representation in that we are not assuming full covariance 
here. Note that in a covariant theory, one usually writes 
down spectral representations in a different form, namely 
S = dM 2 p(M 2 )/{j? L -f x -M 2 ), i.e. with rf-in the 



denominator. This is a special case of Eq. (2.2) with 
p{M 2 lP \) = / °° dM 2 p(M 2 )S(M 2 -M 2 -p\) 



Using thise above ansatz (2J3) for the spectral densi- 
ties, one finds for the self-energy 
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and 



D = x(l - x)f L - xM 2 - (1 - x) (k 2 ± + X 



(2.6) 



Note that vanishes, since 2 jijjji — for j — 

1,2. Self-consistency then requires that 
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j/ L [1 - EzCpI,?!)] +j/ x - [m + £o(pl,pi)] 

(2.7) 

In the above equations we have been sloppy about cut- 
offs in order to keep the equations simple, but this can be 
easily remedied by multiplying each integral by a cutoff 

on the fermion momentum, such as ^A^ — (p — fc)j_ 

In principle, the set of equations [Eqs. (2.4), 



can now be used to determine the spectrum of the model. 
But we are not going to do this here since we are more 
interested in the LF solution to the model. However, we 
would still like to point out that, for large enough g, one 
obtains a self-consistent numerical solution to the Eu- 
clidean version of the model which has a non-vanishing 
scalar piece — even for vanishing current quark mass to, 
i.e. chiral symmetry is spontaneously broken and a dy- 
namical mass is generated for the fermion in this model. 



B. LF Solution 

A typical framework that people use when solving LF 
quantized field theories is discrete light-cone quantization 
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(DLCQ) |l]|. Since it is hard to take full advantage of 
the large Nc limit in DLCQ, we prefer to use a Green's 
function framework based on a 4 component formulation 
of the model. 

In a LF formulation of the model, the fermion prop- 
agator (to distinguish the notation from the one above, 
we denote the fermion propagator by G here) should be 
of the form Q 

G( P n = j + p-G + (2p+p~,f ± )+ 1 -p + G^(2p+p-,f L ) 
+ V±G±(2p + p-,f ± ) + G (2p+p-,f ± ). (2.8) 

Again we can write down spectral representations 

r2 p? F (M\? ± ) 
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where i = +,— ,_L,0. This requires some explanation: 
On the LF, one might be tempted to allow for two terms 
in the spectral decomposition of the term proportional 
to 7 + , namely 



tr(7 _ G) oc / dM 
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(2.11) 



one can cast Eq. (2.10) into the form 
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which is of the form in Eq.(2.9) plus an energy indepen- 
dent term. The presence of such an additional energy 
independent term would spoil the high energy behavior 
of the model p2[ : In a LF Hamiltonian, not all coupling 
constants are arbitrary. In many examples, 3-point cou- 
plings and the 4-point couplings must be related to one 
another so that the high energy behavior of scattering 
via the 4-point interaction and via the iterated 3-point 
interaction cancel [Q. If one does not guarantee such 
a cancellation then the high-energy behavior of the LF 
formulation differs from the high-energy behavior in co- 
variant field theory and in addition one often also gets a 
spectrum that is unbounded from below. In Eq. ( [2.12 ), 
the energy independent constant appears if the coupling 



Note that in a LF formulation, G+ and G_ are not neces- 
sarily the same. 



constants of the "instantaneous fermion exchange" in- 
teraction in the LF Hamiltonian and the boson-fermion 
vertex are not properly balanced. In the following we will 
assume that one has started with an ansatz for the LF 
Hamiltonian with the proper high-energy behavior, i.e. 
we will assu me tha t there is no such energy independent 
piece in Eq. ( |2.12j ). 

The LF analog of the self-energy equation is obt ained 
by starting from an expression similar to Eq.(2.5) and 
integrating over k~ . One obtains 
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and CT is an energy (p~ ) independent counter-term. The 
determination of this counter-term, such that one ob- 
tains a complete equivalence with the Schwinger Dyson 
approach, is in fact the main achievement of this paper. 
First we want to make sure that the counter-term renders 
the self-energy finite. This can be achieved by perform- 
ing a "zero-energy subtraction" with a free propagator, 
analogous to adding self-induced inertias to a LF Hamil- 
tonian, yielding 
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(2.17) 



and where we denoted the finite piece by Am| M (for 
zero-mode), since we suspect that it arises from the dy- 
namics of the zero-modes. Mj 2 is an arbitrary scale pa- 
rameter. We will construct the finite piece (Am 2 ZM ) so 
that there is no dependence on Mq lect in CT in the end. 

At this point, only the infinite part of CT is unique 
P2| , since it is n eeded to cancel the infinity in the k + 
integral in Eq. (2.14), while the finite (w.r.t. the k + 
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integral) piece (i.e. Arn^ M ) seems arbitrary. Q Below 
we will show that it is not arbitrary and only a specific 
choice for Am% M leads to agreement between the SD and 
the LF approach. 

Note that the equation for the self-energy can also be 
written in the form 



sr(p)=s 



dk+ 
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One can prove this by simply comparing expressions! By- 
passing the use of the spectral function greatly simplifies 
the numerical determination of the Green's function in a 
self-consistent procedure. 



T(P) = 



Am% M 



(2.21) 
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(some cutoff, such as a sharp momentum cutoff, is 
implicitly assumed). Having obtained the eigenvalues 
of the DLCQ Hamiltonian, we then determined the 
Gree n's fu nction self-consistently, by iteratively solving 
Eq. ( [2.18| ), using the same cutoffs as in the DLCQ cal- 
culation: the same transverse momentum cutoff and dis- 
crete k + summations instead of the integrals. The result 
was that the invariant mass at the first pole of the self- 
consistently determined Green's function coincides to at 
least 10 significant digits (!) with the invariant mass 
of the physical fermion as determined from the DLCQ 
diagonalization. This result was independent of the cut- 
off — as long as the same cutoff was used in both the 
Green's function and the DLCQ approach. This proves 
that the self-consistent Green's function calculation and 
the DLCQ calculation are in fact completely equivalent 
for our toy model. This is a very useful result, since it al- 
lows us to formally perform the continuum limit (replace 
sums by integrals) — a step that is clearly impossible for 
a DLCQ calculation. 



C. DLCQ solution 



D. Comparing the LF and SD solutions 



There are reasons why one might be sceptical about 
the Green's function approach to the LF formulation 
of the model: First we used a four-component formu- 
lation which resembles a covariant calculation. Fur- 
thermore, we introduced spectral representations for 
the Gre en's functions and assumed certain properties 
[Eq. (|2.9|)1 . Since we were initially also sceptical, we per- 
formed the following calculation: First we formulated the 
above model as a Hamiltonian DLCQ problem jl3| with 
anti-periodic boundary conditions for the fermions and 
periodic boundary conditions for the bosons in the lon- 
gitudinal direction. Zero modes (k + = 0) were omitted. 
This is a standard procedure and we will not give any 
details here. The only nontrivial step was the choice of 
the k ineti c energy for the fermion, which we took, using 
Eq. (pTl6b , to be 



T = E E T (p)(blb p 



(2.20) 



=1,3,. 



with 



Having established the equivalence between the 
Green's function method and the DLCQ approach, we 
can now proceed to compare the Green's function ap- 
proach (in the continuum) with the Schwinger-Dyson ap- 
proach. Motivated by considerations in Ref . M , we make 
the following ansatz for ZM: 
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The motivation for this particular ansatz becomes ob- 
vious one we rewrite the expression for S^ F : For this 
purpose, we first note that 
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Together with the normalization condition 
r o °° dM 2 p FF (M 2 , k 2 ± )=l, this implies 



4 Note that what we called the "finite piece" w.r.t. the k + in- 
tegral is still divergent when one integrates over d 2 k± without 
a cutoff! 
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we considered has an underlying Lagrangian which is par- 
ity invariant, one can use G + = G - for the self-consistent 
solution and still use Eq. ( p.26| ) or Eq. ( [2.27 ) but with 
G + replaced by G_. 
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III. SUMMARY 



(2.24) 



where we used our particular ansatz for Am 2 ZM [Eq. 



(I 



22)]. Thus, for our particular choice for the finite piece 
of the kinetic energy counter term, the expression for 
and T*^ F are almost the same — the only difference be- 
ing the replacement of p+ F with p^_ F and an overall fac- 
tor of p~/p + . Furthermore, the most important result 
of this paper is a direct comparison (take x = k + /p + ) 
shows that the same spectral densities tha t pro vide a self- 
consistent solution to the SD equations (|2.5| ) also yield 
a self-consistent solution to the LF equations, provided 
one chooses 



p I f{M\kl)=p LF {M\ n _ ! 



JjFi-k/tI D\ /A,f2 l2 



p L (M 2 ,k 2 



Po {M\ki) = PQ {M\ki). 



(2.25) 



In particular, the physical masses of all states (in the 
sector with fermion number one) must be the same in 
the SD and the LF framework. 

In the formal considerations above, we found it conve- 
nient to express Am 2 ZM in terms of the spectral density. 
However, this is not really necessary since one can express 
it directly in terms of the Green's function 



Am% M = g 2 p + l dp 



d 2 pi_ 
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Analogously, one can also perform a "zero-energy sub- 
traction" in Eq. (2.18) with the full Green's function, 
i.e. by choosing 



CT = -g 2 J^J^p F G + (2p+p F ,p\ F ) , (2.27) 

with pp = —(A 2 + k 2 _)/2k + . This expression turns out 
to be very useful when constructing the self-consistent 
Green 's fun ction solution. We used both ansatze [Eqs. 
( 2.26 ) and ( 2.27| )1 to determine the physical masses of 
the dressed fermion. In both cases, numerical agreement 
with the solution to the Euclidean SD equations was ob- 
tained. 

Note that, in a canonical LF calculation (e.g. using 
DLCQ) one should avoid expressions involving G+, since 
it is the propagator for the unphysical (" bad" ) compo- 
nent of the fermion field that gets eliminated by solving 
the constraint equation. However, since the model that 



We studied a simple 3+1 dimensional model with 
"QCD-inspired" degrees of freedom which exhibits spon- 
taneous breakdown of chiral symmetry. The methods 
that we used were the Schwinger-Dyson approach, a LF 
Green's function approach and DLCQ. The LF Green's 
function approach was used to "bridge" between the SD 
and DLCQ formulations in the following sense: On the 
one hand, we showed analytically that the LF Green's 
function solution to the model is equivalent to the SD ap- 
proach. On the other hand we verified numerically that 
by discretizing the momentum integrals, that appeared 
in the LF Green's function approach, agreement between 
the LF Green's function approach and DLCQ. Hence we 
have shown that the SD solution and the DLCQ solution 
are equivalent. This remarkable result implies that even 
though the LF calculation was done without explicit zero- 
mode degrees of freedom, its solution contain the same 
physics as the solution to the SD equation — including 
dynamical mass generation for the fermions. 

However, we have also shown that the equivalence be- 
tween the LF approaches and the SD approach only hap- 
pens with a very particular choice for the fermion kinetic 
mass counter-term in the light-front framework. Our cal- 
culation also showed that the current quark mass of the 
SD calculation is to be identified with the "vertex mass" 
in the LF calculation — provided the same cutoffs are 
being used in both calculations. This result makes sense, 
considering that both the current quark mass and the 
LF vertex mass are the only parameters that break chi- 
ral symmetry explicitly. The mass generation for the 
fermion in the chiral limit of the LF calculation occurs 
through the kinetic mass counter-term (which does not 
break chiral symmetry) |hJ .[] Our results contradict Ref. 
Jl4| |, where it has been ad hoc suggested that the renor- 
malized vertex mass remains finite in the chiral limit to 
account for spontaneous breaking of chiral symmetry. 

Our work presents an explicit 3+1 dimensional exam- 
ple that there is no conflict between chiral symmetry 
breaking and trivial LF vacua provided the renormal- 
ization is properly done. In our formal considerations, 
we related the crucial finite piece of the kinetic mass 
counter-term to the spectral density. Several alterna- 
tive determinations (which might be more suitable for a 



5 We should add that the "kinetic mass counter-term" did 
depend on the transverse momentum of the fermion for most 
cutoffs other than a transverse momentum on the fermion. 
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practical calculation) are conceivable: parity invariance 
for physical observables J^] , more input (renormalization 
conditions) such as fitting the fermion or "pion" mass. 

However, one must be careful with this result in the 
following sense: although we have provided an explicit 
example which shows that, even in a 3+1 dimensionsonal 
model with x$B for m — > 0, LF Hamiltonians without 
explicit zero-modes can give the right physics, we are still 
far from understanding whether this is possible in full 
QCD and how complicated the effective LF Hamiltonian 
for full QCD needs to be. More work is necessary to 
answer these questions. 

As an extension of this work we had planned to study 
the pion in the chiral limit of a 1+1 dimensional version 
of this model using the LF framework. [] We were not able 
to derive an analog of the Green's function equations for 
the pion, so we had to resort to a brute force DLCQ cal- 
culation. Numerical convergence, which was acceptable 
for the fermion, was very poor for the pion in the chiral 
limit, and we were thus not able to demonstrate that it 
emerges naturally as a massless particle. Nevertheless, 
we expect that other numerical techniques, which treat 
the end point behavior of the LF wavefunctions more 
carefully than DLCQ, should yield a massless pion for 
this model. 
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Even in 1+1 dimensions, one expects a massless boson in 
the chiral limit because of Nc — > oo. One can show this in the 
SD formalism since the solution for the self-energy equation 
for the fermion also solves the Bethe-Salpeter equation for the 
pseudoscalar bound state with zero mass. 
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